Current-induced torques in the presence of spin-orbit coupling 
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In systems with strong spin-orbit coupling, the relationship between spin-transfer torque and 
the divergence of the spin current is generalized to a relation between spin transfer torques, total 
angular momentum current, and mechanical torques. In ferromagnetic semiconductors, where the 
spin-orbit coupling is large, these considerations modify the behavior of the spin transfer torques. 
One example is a persistent spin transfer torque in a spin valve: the spin transfer torque does not 
decay away from the interface, but approaches a constant value. A second example is a mechanical 
torque at single ferromagnetic-nonmagnetic interface. 
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Introduction — Since the prediction [l|43| of spin trans- 
fer torques in non-collinear ferromagnetic metal circuits, 
they have been the subject of extensive research 0, [5j. 
The possibility of using spin transfer torque to im- 
prove the commercial viability of magnetic random ac- 
cess memory (MRAM) [6], and the rich non-equilibrium 
physics involved establish the topic as one of practical 
and fundamental interest. These torques arise from the 
exchange interaction between non-equilibrium, current- 
carrying electrons and the spin-polarized electrons that 
make up the magnetization. In systems where the spin- 
orbit coupling is weak, the torque on the magnetization 
can be computed from the change in the spins flowing 
through the region containing the magnetization. This 
relation is a consequence of conservation of total spin. 
Here, we consider systems in which the spin-orbit cou- 
pling cannot be neglected (and hence total spin is no 
longer conserved). 

In systems where spin angular momentum is not con- 
served, the relationship between the spin transfer torque 
and the flow of spins needs to be generalized. Conserva- 
tion of total angular momentum implies that mechanical 
torques on the lattice of the material accompany changes 
in the magnetization [3, [8j . This effect has been used for 
decades to measure the ^-factor of metals. More recent 
theoretical [9|, LLQ| and experimental PJJ work considers 
the current-induced mechanical torques present at the in- 
terface of a ferromagnet and non-magnet, similar in spirit 
to the spin transfer torques on the magnetization present 
in spin valves. 

In this article we develop a theory for current-induced 
torques (both spin transfer torques and mechanical 
torques) in systems with strong spin-orbit coupling, and 
apply it to a model of dilute magnetic semiconductors. 
We find that by accounting for the orbital angular mo- 
mentum of the electrons, we can relate the change in 
total angular momentum flow to spin transfer torques 
and mechanical torques. We study two system geome- 
tries where these torques play important roles. The first 
is a spin- valve geometry, which is used to study the fea- 
tures of spin transfer torques in the presence of spin-orbit 
coupling. The second is a single interface between a fer- 
romagnet and non-magnet, which elucidates the physics 
underlying current-induced mechanical torques. 



Formalism — We consider a Hamiltonian consisting of 
a spin-independent kinetic and potential energy Ho = 
~ h 2 ^ + V(r), an exchange splitting A, and an atomic- 
like spin-orbit interaction parameterized by a: 
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where L and s are the electron angular momentum and 
spin operators, respectively [12J. The exchange splitting 
arises from a magnetization M, with magnitude M s . We 
treat the magnetization within mean field theory. 

We consider the torque on the magnetization due to 
electric current flow. The spin transfer torque tstt 
at position r from electronic states with spin density 
s(r) is proportional to the component of spin trans- 
verse to the magnetization [14[: TsxT( r ) = — -j^- = 
^ (M(r) x s(r)). In the absence of spin-orbit coupling, 
this torque can be related to the divergence of a spin 
current, which offers conceptual and computational sim- 
plicity [15j. In the following we analyze how spin-orbit 
coupling changes this simple result. One consequence is 
an expression for the mechanical torque Ti at . 

We develop an expression for tstt by evaluating 
the time-dependence of the electron spin and angular 
momentum densities. To do so, we adopt a Heisen- 

berg picture of time evolution, and evaluate Jf' = 
H : ^(r)0^(r) , where -0(r) is the position operator, 
for the operators O = s, L. This procedure leads to [4(: 
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where Q s (r) = ^(r)v^)s^(r), and the velocity operator 

is given by v = ^ ( V — V ) ; here the arrow superscript 

specifies the direction in which the gradient acts. In ad- 
dition: 
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where Ql(f) = ^(r)| vL + Lv)^(r) (the product 
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FIG. 1: Left and right panels shows GaMnAs band structure 
without and with spin-orbit, respectively (for 72 = 73 = 2.4). 
(arrows indicate spin direction of eigenstates) . The inset 
shows the direction of bulk magnetization, and spin, veloc- 
ity, and k vectors for a single state (in black, red, blue, and 
green). The torque from the misalignment between magne- 
tization and spin equals the torque from the misalignment 
between velocity and k vectors. 



of non-commuting operators L and v is symmetrized). 

We've defined fi at (r) = ^t(r) i^o,L ^(r), which is 

nonzero for a potentials V(r) which break rotational sym- 
metry jl6| . 

We define a total angular momentum J = L + s, a total 



angular momentum current Qj 
Eqs. © and © to obtain: 
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Finally, we take the expectation value of Eqs. (J2]|4]), re- 
placing operators by densities. Eq. ^ is our main formal 
result. When spin-orbit coupling is important, the total 
angular momentum in the conduction electrons couples 
both to the magnetization and the lattice. The coupling 
of electron spin to the lattice requires both spin-orbit cou- 
pling and crystal field potential. The term T\ at changes 
the physical picture of spin transfer torque substantially, 
as is illustrated by considering Eq. (j4]) for a single bulk 
eigenstate: ^ and V • Qj vanish, however tstt and 
Ti a t may both be non-zero, implying a coupling from the 
angular momentum of the lattice to the magnetization. 
This coupling flows from the lattice to the orbital sub- 
system through the crystal field, which then couples to 
the spin through spin-orbit coupling, and finally to the 
magnetization through the exchange interaction. 
Application to DMS - - We apply this general formal- 
ism to a model of a dilute magnetic semiconductor 
(DMS). DMSs are semiconductor host materials which 
become ferromagnetic when doped with magnetic atoms. 



Gai_ x Mn x As is the archetype for these materials, and 
can be described as a system of local moments of Mn 
d-electrons, whose interaction is mediated by holes in 
the semiconductor valence band [17J. The valence states 
are described by the Kohn-Luttinger Hamiltonian H^ L , 
which represents a small-k expansion for a periodic Ho, 
acting in the £ — 1 subspace (describing valence states). 
It is given by: 
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where L are the spin-1 matrices for the p- state orbit als, 
7i, 72, 73 are Luttinger parameters, and k is the Bloch 
wave- vector. Figure Q] shows how the presence of spin- 
orbit coupling affects the band structure. 

For periodic systems the velocity operator can be writ- 
ten as: v = ^^|[-, and spin and angular momentum cur- 
rent densities are again denned as symmetrized products 
of v and L, and v and s. The dynamics of the magne- 
tization occur on a much longer time scale than that of 
the electronic states, so we compute the dynamics from 
a sum over scattering states, for which ^ = ^ = 0. For 
the Luttinger Hamiltonian, the ^-component of fi at is: 
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where the brackets on the second term indicate an an- 
ticommutator. Other components are given by cyclic 
permutation of indices. The first term of Eq. (J6]) can 
be written as v x Ttk = -^ (f x fik). This term can be 
interpreted as a torque on the crystal angular momen- 
tum r x ftk, and results from the misalignment between 
wave vector and velocity. It is generically nonzero for 
any material with a non-spherical Fermi surface. In the 
spherical approximation (72 =73), Eq. (jl|) implies that 
the net flux of total angular momentum into a volume 
is equal to the change of magnetization plus the crystal 
angular momentum inside the volume. 
STT in spin-valves — We first consider a system to study 
the tstt term of Eq. (JH). Figure [2fa) shows the geome- 
try; current flows in the indirection, perpendicular to the 
magnetization of both layers. We focus on the compo- 
nent of torque which is in the plane spanned by the two 
magnetization directions. This in-plane torque is deter- 
mined by the out-of-plane (or z-component) spin density 
[L4J . For the results presented here, we use the param- 
eter values: (71,72,73) = (6.85,2.1,2.9), A = 0.27 eV, 
a = 0.11 eV, Ep = 0.16 eV (E^ is measured from the 
top of the valence band) . The tunnel barrier is described 
by Eqs. (pQ) and ([5|), with A = 0, and with an energy 
offset so that the top of the valence band is 0.1 eV below 
Ey- We calculate the eigenstates numerically and apply 
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FIG. 2: (a)Spin valve geometry: FM layers' magnetization 
points in the x and y (out-of-page) directions, (b) The spin 
transfer torque versus position away from the left normal 
metal-FM interface, which decays to zero in the absence of 
spin-orbit coupling, and does not in its presence, (c) Plot of 
the total spin transfer torque, the net flux of spin current, 
and net flux of total angular momentum current versus FM 
thickness. The linear dependence for large thickness is due to 
a persistent spin transfer torque. 



boundary conditions as described in Ref. [18j . 

Figure [2fb) shows the spin transfer torque density as 
a function of distance away from the interface. We find 
that for a = (no spin-orbit coupling) , the torque decays 
to zero away from the interface, as expected [15]. For 
qj/0, the torque oscillates around a nonzero value, and 
extends into the bulk. Figure [2fc) shows that the total 
spin transfer torque as a function of ferromagnetic (FM ) 
layer thickness Lfm is proportional to thickness for large 
Lfm- This is in contrast to the metallic spin valve, where 
the torque is an interface effect and becomes constant for 
large L FM . 

This persistent spin transfer torque arises because the 
spins of individual eigenstates are not aligned with the 
magnetization (see Fig. [1]) in the presence of spin-orbit 
coupling. The misalignment gives rise to a torque be- 
tween the lattice and the magnetization. In equilib- 
rium, these torques cancel when summed over all occu- 
pied states. However, the presence of a current changes 
the occupation of the bulk states and can give rise to 
a torque yjj [20j in systems without inversion symme- 
try. Inversion symmetry is only very weakly broken in 
bulk GaMnAs, and is not included in the Kohn Luttinger 
Hamiltonian, Eq. (|5]). Here, interfaces between materials 
breaks inversion symmetry. 

The combination of an interface and a current flow 
changes the occupation of the bulk states near the Fermi 
energy (depending on the transmission probabilities of 
individual states across the interface) and induces coher- 
ence between these states. The change in the occupation 
probabilities gives rise to a persistent transverse spin ac- 
cumulation, which only decays through other scattering 
mechanisms not included here (e.g. defect scattering). 



FIG. 3: The total spin transfer torque, the net flux of total 
angular momentum, — 7i a t, and the persistent component of 
the spin transfer torque on a FM layer with Lfm = 30 nm as 
a is increased from to ao = 0.11 eV. 



This spin accumulation gives rise to the persistent spin 
transfer torque. The coherence between the states modi- 
fies the spin accumulation and the torque near the inter- 
face but these corrections decay away from the interface 
due to dephasing. 

Figure [3] shows, as a function of the spin-orbit coupling 
constant a, the values of total spin transfer torque, the 
angular momentum current flux, -Ti at , and the persistent 
contribution to spin transfer torque (for Lfm = 30 nm). 
We determine the persistent contribution from the slope 
of the integrated total versus FM width Lfm at large 
Lfm (see Fig. Et). This procedure neglects the contri- 
butions from coherence near the interface. In this exam- 
ple, the spin transfer torque increases with the addition 
of spin-orbit coupling, largely because of the addition of 
the persistent term. This qualitative behavior depends 
on system parameters: for Ep = 0.34 eV, for example, 
the spin-orbit coupling decreases the total torque. 
Nanomechanical torques in wires — We next consider a 
system which exemplifies that physics of the T\ at term of 
Eq. (jl|): a single interface between GaMnAs and GaAs, 
with the direction of the magnetization parallel to the 
current flow (see Fig. 0^). This is similar to the geom- 
etry considered in previous theoretical and experimental 
work [ol— Till] , The vanishing magnetization in GaAs im- 
plies tstt = 0, so that Ti a t = V ■ Qj, and its total value 
can be deduced from Qj n — Qj ut . We use the same pa- 
rameters as before, except Ep = 0.06 eV, and the top of 
the valence band of both layers coincide. 

Figure 2h shows Ti at in the GaAs layer as a function 
of distance away from the interface (assuming electron 
particle flow from left to right). The total torque (dark 
curve) shows oscillatory decay, while the torque from 
a particular channel (light curve) shows simple oscilla- 
tion. The behavior of the single channel is illustrated in 
Fig. [3)3. We assume specular scattering, so that the in- 
cident state chosen (black circle) transmits into the four 
states of GaAs with equal k x , k y (also shown with black 
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FIG. 4: (a) shows the system geometry. We take electron 
particle flow from left to right, and consider the mechanical 
torque 7i a t in the z direction, (b) shows slices of the Fermi 
surface for the different layers, with (J(k)) superimposed, and 
also shows the J-character of the states specified by the black 
circle. Also shown is the transmission probability for each 
of the states in the GaAs. (c) shows the total mechanical 
torque density in the GaAs as a function of distance from the 
interface (dark curve), and the contribution from the single 
incoming state specified in (b) (dashed curve). 



circles). The character of these states, along with the 
transmission probability, is shown in Fig.HjD. The incom- 
ing state couples most strongly to the state with similar 
J character, but also partially transmits into other states 
with different J character and wave vector k z . These dif- 
ferent scattering channels interfere with each other, lead- 
ing to an oscillatory J(z), with an oscillation period in- 
versely proportional to the splitting of k z wave- vectors of 
the different sheets of the Fermi surface. This splitting is 
from the lattice crystal field and spin-orbit coupling, the 
agents responsible for T\ at . Different channels have dif- 
ferent oscillation periods, so that their total decays away 
from the interface, as happens for spin transfer torques in 
ferromagnets [15]. For the parameters used here, we find 



Q l f z = 1.20^, due to the polarization of the states from 
the magnetization, while Q°j^ = 0.46ft-. Mechanisms 
not considered here, such as spin-flip scattering, ensure 
that Qj? decays to zero away from the interface. 

The mechanical torque is Qy z —Q°J^ = 0.74ft-. For ap- 
propriate experimental conditions, this torque is greater 
than the thermal fluctuations and is a measurable ef- 
fect. We refer the reader to Ref. |9l4llj for details of 
treatment of the torsion dynamics and experimental de- 
tails. The formalism developed here generalizes previous 
work to allow for microscopic evaluation of the electronic 
structure contribution to the current-induced mechanical 
torque. For systems with nonzero magnetization, the mi- 
croscopic form of Ti a t is necessary to determine the parti- 
tioning of total angular momentum flux between torques 
on the magnetization and torques on the lattice. Our the- 
ory neglects other mechanisms of spin relaxation, such as 
disorder-induced spin-flip scattering, so that full calcula- 
tions will require microscopic calculations like these to 
be embedded in diffusive transport calculations. 

Conclusion — We have shown how atomic-like spin-orbit 
coupling affects current-induced torques: both the spin 
transfer torque on the magnetization and the mechanical 
torque on the lattice. In GaMnAs spin valves, we find 
a contribution to the spin transfer torque that persists 
throughout the bulk. This result may explain experi- 
ments which find critical currents which are up to an 
order of magnitude smaller than the value expected from 
a simple accounting of the net spin current flux |2l|, [22j . 
For a single interface between GaMnAs and GaAs, we 
microscopically compute the mechanical torque due to 
scattering from the interface. These results highlight im- 
portant, qualitatively different physics at play when spin- 
orbit coupling is strong. 

The authors acknowledge helpful conversations with A. 
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